Critical behavior in reaction-diffusion systems exhibiting absorbing phase transition 
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Phase transitions of reaction-diffusion systems with site occupation restriction and with particle 
creation that requires n > 1 parents and where explicit diffusion of single particles {A) exists are 
reviewed. Arguments based on mean-field approximation and simulations are given which support 
novel kind of non-equilibrium criticality. These are in contradiction with the implications of a 
suggested phenomenological, multiplicative noise Langevin equation approach and with some of 
recent numerical analysis. Simulation results for the one and two dimensional binary spreading 
2 A — > 4j4, 4j4 — > 2 A model display a new type of mean-field criticality characterized by a = 1/3 
and P = 1/2 critical exponents suggested in cond- mat/0210615. 
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I. INTRODUCTION 

The classification of universality classes of second or- 
der phase transitions is still one of the most important 
uncompleted task of statistical physics. Recently phase 
transitions of genuine nonequilibrium systems have inten- 
sively been investigated among reaction-diffusion (RD) 
type of models exhibiting absorbing states [1-3]. There 
has been a hope that in such homogeneous systems sym- 
metries and spatial dimensions are the most significant 
factors like in equilibrium ones, but gradually it turned 
out that there may be more relevant constituents as well. 
The most well known example is the parity conserving 
class (PC), which differs from the robust universality 
class of directed percolation (DP). The DP hypothesis 
stated by Janssen and Grassberger [4,5], according to 
which in one component systems exhibiting continuous 
phase transitions to single absorbing state (without extra 
symmetry and inhomogeneity or disorder) short ranged 
interactions can generate DP class transition only. How- 
ever parity conservation itself proved to be an insufficient 
condition in many cases [6-9] and rather an underlying 
A 3A, 2A (I) (BARW2) process [10] of particles 
and the Z2 symmetry of absorbing states is necessary 
for this class [11]. On the other hand parity conserva- 
tion in N-component branching and annihilating random 
walk (N-BARW) systems [10], or by triplet production 
models [12] was found to be responsible for novel classes 
again. In one dimensional, multi-component reaction- 
diffusion systems site restriction turned out to be a rel- 
evant, new factor [13,14]. Global conservation laws by 
directed percolation and lattice gas models were shown 
to be irrelevant [15-18], while systems with multiple ab- 
sorbing states [19] or with multi-components also exhibit 
DP class scaling behavior [20] . 

An other important puzzle was being investigated in- 
tensively during the past three years that emerges at 
phase transitions of binary production reaction-diffusion 
systems [8,9,20-29] (PCPD). In these systems particle 
production by pairs competes with pair annihilation and 
single particle diffusion. If production wins steady states 
with finite particle density appear in (site restricted), 
while in unrestricted (bosonic) models the density di- 



verges. By lowering the production/annihilation rate a 
doublet of absorbing states without symmetries emerges. 
One of such states is completely empty, the other pos- 
sesses a single wandering particle. In case of site re- 
stricted systems the transition to absorbing state is con- 
tinuous. It is important to note that these models do not 
break the DP hypothesis, because they exhibit multiple 
absorbing states which are not frozen, lonely particle(s) 
may diffuse in them. However no corresponding symme- 
try or conservation law has been found yet. Non-DP type 
of phase transition in a binary production system was al- 
ready mentioned in the early work of Grassberger [30]. 
A corresponding bosonic field theoretical model the an- 
nihilation fission (AF) process was introduced and stud- 
ied by Howard and Tauber [21]. These authors claim 
a non-DP type of transition in AF, because the action 
does not contain linear mass term and the theory is non- 
renormalizible perturbatively unlike the Reggeon field 
theory of DP. In field theories of models exhibiting DP 
class transition the canonical A ^ A 2A reactions 
are generated by the renormalization transformation un- 
like here. Further facts opposing DP criticality are the 
set of different mean-field exponents and the different 
upper critical dimension of binary production PCPD like 
models (4 = 2 vs. 4 = 4) [21,9]. 

Forthcoming numerical studies reported somewhat dif- 
ferent critical exponents, but there has been a consensus 
for about two years that this model should possess novel, 
non-DP type of transition. The first density matrix study 
by Carlon et al. [22] did not support a DP transition, but 
reported exponents near to those of the PC class. Since 
the PCPD does not conserve the particle number modulo 
2, neither exhibits Z2 symmetric absorbing states the PC 
criticality was unfavored. Simulation studies by Hinrich- 
sen [23] and Odor [24] and coherent anomaly calculations 
by Odor [24,26] resulted in novel kind of critical behavior, 
although there was an uncertainty in the precise values of 
critical exponents. Exponent estimates showed diffusion 
(D) dependence that was under-pinned by pair mean- 
field results [24], possessing two distinct classes as the 
function of D. Recently Park et al. reported well defined 
set of critical exponents in different versions of binary 
production PCPD-hke processes [31] . However these sim- 
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ulations were done at a fixed, high diffusion/annihilation 
rate and agree with Odor's corresponding results [32]. 
Kockelkoren and Chate on the other hand claim an other 
set of critical exponents [33] that agrees with Odor's low 
diffusion/annihilation data. 

The PCPD model can be mapped onto a two- 
component model [25] in which pairs are identified as 
a particle species following DP process and single par- 
ticles as an other, coupled species following annihilat- 
ing random walk. Simulations of such a two-component 
system at D = 0.5 showed a continuous phase transi- 
tion with exponents agreeing with those of the PCPD 
for high diffusions. This model is similar to an other one 
[34], which exhibits global particle number conservation 
as well. Field theory [34] and simulations [35,36] for the 
latter model reported two different universality classes as 
the function of D. It would be interesting to see if this 
conservation law is relevant or not like in case of the DP 
[18]. 

Interestingly, higher level cluster mean-field approxi- 
mations result in a single class behavior by varying D and 
it turned out that by assuming logarithmic corrections 
the single class scenario can be supported by simulations 
too [32]. The origin of such logarithmic corrections may 
be a marginal perturbation between pairs and single par- 
ticles in a coupled system description. A filed theoretical 
explanation would be necessary. 

Two more recent studies [37,38] reported non- 
universality in the dynamical behavior of the PCPD. 
While in the former one Dickman and Menezes explored 
different sectors (a reactive and a diffusive one) in the 
time evolution and gave non-DP exponent estimates, in 
the latter one Hinrichsen set afloat a speculative con- 
jecture that the ultimate long time behavior might be 
characterized by DP scaling behavior. In a forthcoming 
preprint [39] Hinrichsen provided a discussion about the 
possibility of the DP transition based on a series of as- 
sumptions. His starting point is a Langcvin equation that 
is mapped onto a wetting process by Cole-Hopf transfor- 
mation. By analyzing this process within a certain po- 
tential he gave arguments for a DP transition. While this 
Langevin equation with real noise is valid for the bosonic 
version of PCPD at and above the critical point, its usage 
in case of site occupancy restricted models is hypothetic, 
the noise can be complex at the transition point and may 
even change sign by the transformation. Furthermore the 
diffusive field of solitary particles is neglected. 

In a very recent preprint [40] Barkema and Carlon 
continue this line and show that some simulation and 
density matrix renormalization results may also be in- 
terpreted as a signal of a phase transition belonging to 
the DP class. By assuming correction to scaling expo- 
nents that are equal to DP exponents and relevant up 
to quadratic or 3-rd order in the asymptotic limit they 
fitted their numerical results in case of two independent 
exponents. The extrapolations resulted is close to DP 
values for D = 0.5. However for smaller D-s and by 



surface critical exponents this technique gave exponent 
estimates which are out of the error margin of DP. 

An other novel class that may appear in triplet pro- 
duction systems was proposed in [33,41] (TCPD). This 
reaction-diffusion model differs from the PCPD that for 
a new particle generation at least three particles have to 
meet. For such generalizations Park et al. proposed a 
phenomenological Langevin equation that exhibits real, 
multiplicative noise [41]. By simple power counting they 
found that the triplet model exhibits distinct mean-field 
exponents and upper critical dimension 4/3 < dc < 
8/3. The simulations in Id [41] indeed showed non- 
trivial critical exponents, which do not seem to corre- 
spond to any known universality classes. Kockelkoren 
and Chate reported similar results in stochastic cellular 
automata (SCA) versions of general nA {n + k)A, 
mA — > {m — l)A type of models [33], where multiple 
particle creation on a given site is suppressed by an ex- 
ponentially decreasing creation probability {p^^^) of the 
particle number. They claim that their simulation results 
in Id are in agreement with the fully occupation number 
restriction counterparts and set up a general table of uni- 
versality classes, where as the function of n and m only 4 
non-mean-field classes exist, namely the DP class, the PC 
class, the PCPD and TCPD classes. However more ex- 
tensive simulations of 1 and 2 dimensional site restricted 
lattice models [12] do not support some of these results 
in case of different triplet and quadruplet models. In 
the 3^ 4A 3^ ^ triplet model Id numerical data 
can be interpreted as mean-field behavior with logarith- 
mic corrections and in two dimensions clear mean-field 
exponents appear, hence the upper critical dimension is 
dc = 1, which contradicts the Langevin equation predic- 
tion. Surprisingly other non-trivial critical behavior were 
also detected in the 'SA 5A 3A A parity conserving 
triplet model and in some quadruplet models [12]. The 
cause of differences between the results of these studies 
is subject of further investigations. Again proper field 
theoretical treatment would be important. 

The classification of universality classes of nonequilib- 
rium systems by the exponent /x of a multiplicative noise 
in the Langevin equation was suggested some time ago 
by Grinstein et al. [42] . However it turned out that there 
may not be corresponding particle systems to real multi- 
plicative noise cases [21] and an imaginary part appears 
as well if one derives the Langevin equation of a RD 
system starting from the master equation in a proper 
way. Furthermore for higher-order processes the emerg- 
ing nonlinearities in the master equation action do not 
allow a rewriting in terms of Langevin-type stochastic 
equations of motion, hence for high-order processes like 
those of the TCPD a Langevin representation may not 
exist. 

This situation resembles to some extent to a decade 
long debate over the critical phase transition of driven 
diffusive systems [43 45]. The latest papers in this topic 
suggest that the phenomenological Langevin equation 
originally set up for such systems do not correspond ex- 
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actly to the lattice models investigated. Simulations for 
different lattice models show, that instead of an exter- 
nal current the anisotropy is the real cause of the critical 
behavior observed in simulations [46,47]. 

II. MEAN-FIELD CLASSES 

In this section I show that mean-field classes of site re- 
stricted lattice models with general microscopic processes 
of the form 

nA^{n + k)A, mA^ {m-l)A, (1) 

with n > 1, TO > 1, A: > 0, Z > and to — / > are 
different from those of the DP and PC processes backing 
numerical results which claim novel type of criticality 
below dc- The mean- field equation that can be set up 
for the lattice version of these processes (with creation 
probability a and annihilation or coagulation probability 
A = 1 — ct) is 

^=akap"{l-p)''-al{l-a)p^, (2) 

where p denotes the site occupancy probability and o is a 
dimension dependent coordination number. Each empty 
site has a probability (l-p) in mean- field approximation, 
hence the need for k empty sites at a creation brings in 
a (1 — p)*^ probability factor. By expanding (1 — p)'^ and 
keeping the lowest order contribution one can see that 
for site restricted lattice systems a p"+^-th order term 
appears with negative coefficient that regulates eq.(2) 
in the active phase. In the inactive phase one expects 
a dynamical behavior dominated by the mA — > pro- 
cess, for which the particle density decay law is known 
p{t) oc [10]. The steady state solutions were 

determined in [12] analytically and one can distinguish 
three different situations at the phase transition: (a) 
n = TO, (b) n > TO and (c) n < to . 

A. The n = m symmetric case 

As discussed in [12] the leading order singularities of 
steady state solution can be obtained. By approaching 
the the critical point CTc = in the active phase the 
steady state density vanishes continuously as 

poc [ct-CTcI'^ , (3) 

with the order parameter exponent exponent = 1. 

At the critical point the density decays with a power-law 

p (X t-"'"" , (4) 

with aMF = (3^^ /vy^^ = 1/n, hence = n, provid- 

ing a series of mean-field universality classes for n > 1 
(besides DP an PC where Vy^^ = 1) and backing the 



results, which claim novel type of non-trivial transitions 
below the critical dimension. Unfortunately determining 
the the value of dc is a non-trivial task without a proper 
Langevin equation. These scaling exponents can be ob- 
tained from bosonic, coarse grained formulation too [41], 
where a p""'"^-th order term, with negative coefficient had 
to be added by hand to suppress multiple site occupancy. 
It is known however that hard-core particle blocking may- 
result in relevant perturbation in d = 1 dimension [14], 
so for cases where the upper critical dimension is dc > 1 
the site restricted, N > 1 cluster mean-field approxima- 
tion that takes into account diffusion would be a more 
adequate description of the model (see [48]). 



B. The n > m case 



In this case the mean-field sohition provides first order 
transition (see [12]), hence it does not imply anything 
with respect to possible classes for models below the crit- 
ical dimension {d < dc). Note however, that by higher 
order cluster mean-field approximations, where the diffu- 
sion plays a role the transition may turn into continuous 
one (see for example [49-51]). The simulations by Kock- 
elkoren and Chate report DP class transition for such 
models in one dimension [33]. 



C. The n < m case 



In this case the critical point is at zero branching rate 
(Tc = 0, where the density decays with a^^ = 1/(to — 1) 
as in case of the n = 1 branching and m = I an- 
nihilating models showed by Cardy and Tauber [10] 
(BkARW classes). However the steady state solution 
for particle production with n > 1 parents gives differ- 
ent f] exponents than those of BkARW classes, namely 
pMF _ _ defining a whole new series of mean- 

field classes [12], for a simplicity I shall call them PkARW 
classes. It is important to note, that one has not found 
a corresponding symmetry or conservation law to these 
classes. These mean-field classes imply novel kind of crit- 
ical behavior for d < dc- For n = 2, m = 3 the mean- 
field exponents are /3 = 1 and a = 1/2 agreeing with 
the mean-field exponents of the PCPD class. Indeed for 
n ~ 2, m = 3 [33] reports PCPD class dynamical criti- 
cality. This supports the expectation that non-mean-field 
classes follow the distinctions observed in the correspond- 
ing mean-field classes. To go further in sections III and 
IV I investigate the phase transition of the simplest un- 
explored PkARW classes, in the n = 2, to = 4 model. 
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D. The role of k and I 



In the mean-field approximation k and I do not af- 
fect the universal properties, however simulations in one 
dimension showed [12] that in case of the to = ri = 
k = I — Z model the critical point was shifted to zero 
branching rate and a BkARW class transition emerged 
there contrary to what was expected for n — m. For a 
stochastic cellular automaton version of these reactions 
[33] reported a non-trivial critical transition. In general 
one may expect such effects for large k and I values, for 
which N-cluster mean-field approximation - that takes 
into account diffusion - would give a better description. 
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FIG. 1. in the two dimensional 2A — > 4A, 

AA 2A model at Gc = 0. The dashed line shows a lin- 
ear fitting for t > 10^ MCS resulting in a = 0.334(1). 



III. SIMULATIONS OF THE 2A 4A, 4A ^ 2A 
MODEL IN TWO DIMENSIONS 



In the II B section I introduced PkARW mean-field 
classes for n < m. Here I explore the phase transition 
in the simplest model from this class, in the 2A AA, 
4A 2A model with D = 0.5 difltusion rate. Two di- 
mensional simulations were performed on i = 1000 lin- 
ear sized lattices with periodic boundary conditions. One 
Monte Carlo step (MCS) — corresponding to dt — 1/P 
(where P is the number of particles) — is built up 
from the following processes. A particle and a number 
X e (0, 1) are selected randomly; if x < D = 0.5 a site 
exchange is attempted with one of the randomly selected 
empty nearest neighbors (nn); ii x > D = 0.5 two parti- 
cles are created with probability a at randomly selected 
empty nn sites provided the number of nn particles was 
greater than or equal 2; or if x > 0.5 two particles are 
removed with probability 1 — a. The simulations were 
started from fully occupied lattices and the particle den- 
sity {p{t)) decay was followed up to 4 x 10^ MCS. 

First the critical point was located by measuring the 
dynamic behavior of p{t). It turned out that the tran- 
sition is at zero branching rate (cc = 0). The density 
decay was analyzed by the local slopes defined as 



Density decays for several as in the active phase 
(0.0002 < (J < 0.05) were followed on logarithmic time 
scales and averaging was done over ~ 100 independent 
runs in a time window, which exceeds the level-off time by 
a decade. The steady state density in the active phase 
near the critical phase transition point is expected to 
scale as 

p(oo, fj) oc |cr - acf ■ (6) 

0.55 I ' , ' , ' 1 
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FIG. 2. /9e// as the function of a^^ in the two dimensional 
2 A -> 4A, AA 2 A model (bullets). The dashed line shows 
a fitting of the form (8). 



-\n[p{t)/p{t/7n)] 



(5) 



where I used to = 4. As Fig.l shows the local slopes 
curve for i > 10^ MCS extrapolates to the mean-field 
value a = 0.334(1). This value agrees with the mean- 
field value a*^^ = 1/3. 



Using the local slopes method one can get a precise 
estimate for (3 as well as for the corrections to scaling 



lnp((X), aj) - lnp(oo, (Tj^i) 
In(cri) - ln((Ti_i) 



(7) 



As one can see on Fig. 2 the effective exponent clearly 
tends to the expected mean-field value (3 = 0.5 as cr ^ 0. 
Assuming a correction to scaling of the form 
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P,ff=P-at-^ (8) 

non-linear fitting results in 5 = 0.44(1) correction to scal- 
ing exponent. 

Besides these scaling correction assumptions I also 
tried to apply different, lowest order logarithmic cor- 
rections to the data, but these fittings gave exponents 
slightly away from mean-field values and the correspond- 
ing coefficients proved to be very small, therefore I con- 
cluded that dc < 2. In the next section I do the same 
analysis in d = 1. 



IV. SIMULATIONS OF THE 2A 4A, 4A ^ 2A 
MODEL IN ONE DIMENSIONS 

The simulations in one dimension were carried out on 
L — 20000 sized systems with periodic boundary condi- 
tions. The initial states were again fully occupied lattices, 
and the density of particles is followed up to 4 x 10^ MCS. 
An elementary MCS consists of the following processes: 

(a) yl0 ^ 9 A with probabihty D, 

(b) AAAA with probability (1 - a)(l - D), 

(c) AAM AAAA or 00^^ AAAA with probabihty 

The critical point was located at (Tc = again. As one 
can see on Figure 3 there is a crossover of the local slopes 
for t > 5 X 10^ MCS and a linear extrapolation for this re- 
gion results in a = 0.329(5) agreeing with the mean-field 
value. 

0.2 I ' 1 ' : ' : ' 1 




0.4 I ' < ' < ' < ' 1 

0.000 0.010 0,^0 0.030 0.040 

FIG. 3. Qe//(t~^^^) in the one dimensional 2A ~* 4A, 
4 A 2 A model at (Jc = 0. The dashed line shows a lin- 
ear fitting for t > 5 X 10^ MCS resulting in a = 0.32(2). 

The steady state data were analyzed in the active re- 
gion for 0.0003 < a < 0.5 as in two dimensions, by the 
local slopes method (eq. 7) and by assuming correction 
to scaling of the form (8). This resulted in 6' = 0.332 



correction to scaling exponent. The local slopes plot- 
ted as the function shown on Fig. 4 extrapolates to 
(3 = 0.49(1) in agreement with the mean- field exponent. 
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FIG. 4. Peff as the function of a"* in the one dimensional 
2 A 4A, 4A 2 A model (bullets). The dashed line shows 
a fitting of the form (8). 



V. CONCLUSIONS 

In this paper I reviewed and discussed the state of the 
art of the phase transitions of reaction-diffusion systems 
exhibiting explicit diffusion and production by n > 1 par- 
ents. Arguments are given against DP criticality that has 
recently been suggested in some papers. These are sup- 
ported by a series of mean-field classes that can be classi- 
fied by the existence of a n = m symmetry in the system 
and by the n and m values. Especially the need for a 
proper field theoretical treatment is emphasized. The 
upper critical dimension in these models is not known. 

I determined by simulations the a and f3 exponents of 
the 2A — > AA, AA 2A model in one an two dimen- 
sions. These results indicate that for this binary produc- 
tion system the upper critical dimension is dc < 1. This 
model conserves the parity of particles still its transition 
does not belong to the PC class. Note that while the den- 
sity decay results in one dimension are in agreement with 
that of Kockelkoren and Chate-s simulations [33] , the off- 
critical order parameter exponent is f3^'^^ — l/(m — n) 
which shows that there are more classes exist at zero 
branching rate besides the BkARW universality classes. 
It is still an open question if there is any variant of 
PkARW models that exhibits non-mean-field criticality 
in physical dimensions {d > 1). 
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